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Abstract 

The chiral two-pion exchange component of the long-range pp interaction 

is studied in an energy-dependent partial-wave analysis. We demonstrate its 

presence and importance, and determine the chiral parameters q (i = 1,3,4). 

The values agree well with those obtained from pion-nucleon amplitudes. 
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The longest-range part of the strong nucleon-nucleon (NN) interaction is the well- 
established one-pion exchange (OPE) force Next in range is the two-pion exchange 
(TPE) force, the formulation of which has been a long-standing problem both in field 
theory [f|,|5| and in dispersion theory 0. In recent years, it has been argued that the key to 
the solution is the chiral symmetry of QCD [7J, and that the long-range parts of the TPE 
potential can be derived model-independently by a systematic expansion of the effective 
chiral Lagrangian ||. In this Letter, we will study this long-range chiral TPE force in the 
proton-proton (pp) interaction and show unambiguously its presence and its importance. 

In the energy-dependent Nijmegen partial-wave analyses (PWA's) of the NN and NN 
scattering data 0-0, the long-range forces are taken into account exactly and the short- 
range forces are parametrized analytically. The partial-wave scattering amplitudes are ana- 
lytic functions of the energy. The nearby left-hand singularities in the complex-energy plane 
are due to the long-range forces; these cause the rapid energy dependence of the physical 
NN scattering amplitudes. The shorter-range forces are responsible for the far-away sin- 
gularities, which give in the physical region only slow energy variations of the amplitudes. 
This method of PWA can serve as a sensitive tool to investigate precisely these long-range 
interactions. It has been used successfully in studies of electromagnetic interactions ]T3|] and 
of the OPE potential Here this tool will again be employed, now to study the 

long-range chiral TPE component of the pp force. 

The methods of the Nijmegen PWA's are described in detail in Refs. f|-[l2|- The 
long-range potentials, including the full electromagnetic interaction (relativistic Coulomb, 
magnetic-moment interaction, and vacuum polarization) and the longest-range strong inter- 
actions are used in the relativistic Schrodinger equation which is solved with a boundary 
condition (BC) at some r = b. This BC is parametrized as an analytic function of energy 
for the various partial waves. The BC parameters, representing short-range physics, and 
the free parameters in the long-range forces (e.g. the pion-nucleon coupling constant) are 
determined from a fit to the data. In the "standard" Nijmegen PWA's of Refs. |rO| , |TTf the 
boundary is put at b = 1.4 fm, and the long-range strong potential outside of 1.4 fm is 



taken as the OPE potential supplemented by the non-OPE forces of the Nijmegen soft-core 
potential Nijm78 [|17]]. These heavy-boson exchanges were included because OPE alone did 
not allow for an optimal description of the data. In this standard pp PWA, we obtain 
with 19 BC parameters xLn = 1968.7 and f^ = 0.0756(4), where the error is statistical, 
on the Nijmegen 1998 pp database below 350 MeV, in which 1951 pp scattering data are 
included ||18|| . This result will serve here as a benchmark. 

Let us demonstrate our method first with some parts of the electromagnetic interaction. 
When one omits in the standard 1998 pp PWA the magnetic-moment interaction, both from 
the potential and in constructing the scattering amplitude, the xLm increases by 390.0 to 
Xmin = 2358.7. This is therefore a 19.7 standard deviation (s.d.) effect. Omitting vacuum 
polarization leads to xLm = 2181.3, i.e., a rise in xLm °f 212.6, which corresponds to 14.6 
s.d. These numbers demonstrate that one can use this method of energy-dependent PWA to 
show the presence and the importance of these specific well-known parts of the long-range 
pp interaction. 

A very important part of the energy dependence of the AW phase shifts comes from OPE. 
In the Nijmegen energy-dependent PWA's the different pion-nucleon coupling constants 
could be determined accurately and reliably [Pt |l4"| , |T5T| . In Ref . [|T5| , we recommended for the 
charge-independent coupling constant the value f% NlT = 0.0750(9), where the error includes 
statistical as well as systematic effects. As a systematic check, the masses of the exchanged 
pions were determined, with excellent results: m n o = 135.6(1.0) MeV and m n + = 139.6(1.3) 
MeV. In this way, the presence of OPE in the AW force was shown with an enormous 
statistical significance. A more subtle effect is the energy dependence of the OPE potential 
due to the minimal-relativity factor M/E, where M is the proton mass and E the proton 
center-of-mass energy. Omitting this factor from the OPE potential results in Xmm = 1977.2. 
This is a rise of 8.5 in Xmin? or an almost 3 s.d. effect. Recently, also the electromagnetic 
corrections to the OPE potential in np scattering were investigated [T^ . 

The starting point to derive the OPE and TPE potentials is the effective chiral La- 
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jrangian, the leading-order of which is the nonlinear Weinberg model [ITS 



jgf(o) = _]v 



N 



with the chiral-covariant derivatives 



(1) 



^N= (dP + yCoT-TtxDI^N . (2) 

Here, D = 1 + vr 2 /i^, fiU = 1-2573 is the Gamow- Teller coupling, and = 185 MeV is 
the pion decay constant; chiral symmetry fixes cq = 1. Eq. (P implies that the planar- 
and crossed-box TPE diagrams should be calculated with the pseudovector (PV) NNtt 
Lagrangian. We use the physical NNn coupling constant /, i.e., we trade in the Goldberger- 
Treiman value g&j F n for y/Anf /m s ; the scaling mass m s serves to make / dimensionless and 
is conventionally chosen to be numerically equal to the charged-pion mass, m s = m n +. In 
addition to the PV NNtt interaction, Eq. (JI|) contains the Weinberg- Tomozawa (WT) NN2tt 
seagull interaction [GO], resulting in triangle and football TPE diagrams. 

In order to derive the TPE potential in subleading order, three more NN2tt interactions 
are required M, viz. 



JSf« 



-N 



8c 1 D- 1 mir/Ft + 4c a D lt .D'- 



+ 2c 4 t ■ x D l 



N, 



(3) 



leading to additional triangle diagrams. The values of the chiral parameters ("low-energy 
constants") q (i = 1, 3, 4) of order (1/M) are not fixed by chiral symmetry; the q's represent 
"integrated-out" hadrons, such as the heavier mesons like the e and g, and the N- and 
A-isobars. The definition Eq. (^) of these q's [H] agrees with the convention used in heavy- 
baryon %PT ; an additional C2-term does not contribute to the NN force in this order. 

The ci-term violates chiral symmetry explicitly. A systematic expansion of Eqs. ([!]) and ([|) 
to order (1/M) gives the relevant part of the chiral Lagrangian . 



4 



The OPE and TPE potentials derived from this Lagrangian contain central, spin-spin, 
tensor, and spin-orbit terms, viz. 

V = V c + V S (t 1 -<t 2 + V t S 12 + V S oL- S , (4) 

where Si 2 = 3 <j\ ■ r cr 2 ■ r — cti ■ cr 2 . With i = C, S, T, SO, £ = m n / m s , and x = m w r, we 
can write 

Vt(r) = f 2n e n (M/E) [ Vi (x) + w i (x)r v r 2 \m n , (5) 

where n = 1 for OPE and n = 2 for TPE. 

The long-range OPE potential contains an isovector spin-spin part ws and an isovector 
tensor part %, 

ws(x) = e~ x /3x , 

w T {x) = (l + x + x 2 /3) e~ x /x 3 . (6) 

For the pp case, the neutral-pion mass m^o is used in OPE. The coupling f£ = f^ pn0 is a 
free parameter. 

For TPE, the dimensionless isoscalar functions i>, are written as the sum of the leading- 
order terms v it i and the subleading-order terms v^ 2 , 

Vi(x) = (2/tt) v^x) + (m n /M) v i;2 {x) , (7) 

and similarly for the isovector functions In the TPE potential, we use the average 
pion mass m n = 138.04 MeV and the fixed charge-independent coupling constant is f 2 = 
Inn-k = 0.0750. Care must be taken to obtain the appropriate form for the use of Eq. (|5|) in 
the relativistic Schrodinger equation. Other forms of the OPE potential or other two-body 
equations will, in general, give different TPE potentials P,p5]. 

The leading-order static potential TPE(l.o.) contains isoscalar spin-spin and tensor 
terms, vs,i and vt,\ respectively, and an isovector central component u>c,i- The long-range 
parts are 
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v s ,i(x) = 12K (2x)/x 3 + (12 + 8a; 2 )Ai(2a;) / ' x 4 , 

v T ,i(x) = -12K (2x)/x 3 - (15 + 4:r 2 )A 1 (2a;)/2 4 , (8) 
w c ,i{x) = (c 2 + 10c - 23 - Ax 2 ) K (2x)/x 3 + 



(c 2 + 10c - 23 + (45 - 12)x 2 ) K x {2x)/ 



x 4 



where the modified Bessel functions have asymptotic behavior K n (2x) ~ a/vt/4x e 2x . This 
TPE(l.o.) is the "TMO" potential [p6 |, supplemented by the diagrams with the WT seag- 



ulls |2^j2^|. In the WT terms we extracted, for ease of presentation, an overall factor f 4 , cf. 
Eq. (H), and defined c = c /g\, where §a = F v \fA~Kf/m s . 

The subleading-order potential TPE(s.o.) contains nonstatic terms from Eq. ([!]) and the 
leading-order terms from Eq. (Q). The long-range parts read 

M*) = E°=i S e _2 7^ , (9) 

and similarly for 10^2, with the coefficients a p as collected in Table [I]. Also here a factor f 4 
was extracted and the result was rewritten in terms of Co and q = CiM/g\. Our results for 
TPE(s.o.) agree with Ref. [|| . 



Remarkably, a large part of the correct TPE potential was already obtained by Sugawara 
and Okubo in "pre-chiral days," by using PV coupling and two phenomenological NN2n 
interactions: the WT term of Eq. (0) and the ci-part of Eq. @. They also pointed out that 
PV coupling gives a rather strong attractive isoscalar spin-orbit force in subleading order. 
However, the important additionally required chiral C3- and C4-terms were missing; these 
were for the AW case first given in Ref. ||. 

We now come to the results of the TPE studies, in which we again use the 1998 database 
below 350 MeV, with 1951 data |18[| . The main results of the various PWA's are summarized 
in Table ||. We start conservatively with the boundary at b = 1.8 fm, since beyond 1.8 fm 
only OPE and TPE are expected to contribute significantly. When only OPE is included 
as strong force, xLm = 1956.6 is reached at the cost of 29 BC parameters. We want to 
investigate if the fit can be even further improved when TPE is added. When only the 



TPE(l.o.) potential of Eq. (|j) is used, we obtain xLin = 1965.9 with 26 BC parameters. But 
we can do better. The complete TPE potential, xTPE = TPE(l.o.) + TPE(s.o.), contains 
three a priori unknown constants: the chiral parameters q (i = 1,3,4) from Eq. (|3]). In 
the fits we obtain c\ = — 4.4(3. 4)/GeV. The values of c\ and C3, appearing both only in the 
isoscalar central potential, cf. Table Q, are strongly correlated. The correlations between the 
parameters can be summarized concisely by: 

c 3 = [-5.08 - 0.62(ci + 0.76) + 40(/ 2 - 0.0755)] /GeV, 
c 4 = [+4.70 + 0.01(ci + 0.76) + 250(/ 2 - 0.0755)] /GeV. 



In order to determine reliable values for C3 and C4, we use the theoretical estimate |22| for 
Ci obtained from the scalar form factor a(t) of the proton at t — 0, viz. 

Cl = -[a(0)/4m 2 + 9A 2 /16^] ; (10) 

cr(0) is the pion-nucleon sigma term, the value of which is uncertain. We take here the 



plausible "low" value cx(0) = 35(5) MeV | 31j| , which is supported by the recent tcN PWA of 
Ref. [j32] . This gives 



Cl = -[0.46(7) + 0.30]/GeV = -0.76(7)/GeV ; (11) 

the error here is theoretical. Our determination of c\ is consistent with this value. Fixing 
Cl = -0.76/GeV, we find, with 22 BC parameters, xLn = 1937.8 and / 2 = 0.0755(7); the 
resulting values for C3 and C4 are 

c 3 = -5.08(28)/GeV , c 4 = +4.70(70)/GeV , (12) 

where the errors are statistical. The improvement over only OPE is reflected, even beyond 
1.8 fm, in the 18.8 lower Xmm an d i n the 7 fewer BC parameters required. 

The result found for / 2 is in very good agreement with the value 0.0756(4) determined 
in the standard 1998 pp PWA. Our values for the q's can be compared to the determination 



from the irN scattering amplitudes in Ref. |$3|. Here, C\ = — 0.93(9)/GeV was obtained 



using Eq. (0), but with a(0) = 45(8) MeV, along with c 3 = -5.29(25)/GeV and c 4 = 
+3.63(10)/GeV. In view of the uncertainties in the 7rN amplitudes [Q, the good agreement 
is a significant success. It underlines, for the first time quantitatively, that the long-range 
NN and the low-energy ttN interactions are governed by the same chiral Lagrangian. 

In previous studies of the OPE potential, a good systematic check has been the deter- 
mination of the masses of the exchanged pions. In order to check explicitly that we are now 
actually looking at the TPE interaction, we determine the range. This is done by adding 
the pion mass m n in the potential ;^TPE as another free parameter. We first fix the pion 
coupling in OPE at f£ = 0.0755 and the q's to their central values given in Eqs. (|TTD and 
(|I~2|). Then we fit an overall scale factor A for the potential %TPE, the pion mass m n , and the 
BC parameters. The results are: A = 0.82(16) and = 125(10) MeV. Alternatively, we fix 
c\ and fit m n together with C3, C4, and the BC parameters. This results in m n = 128(9) 
MeV, again in good agreement with the average pion mass = 138.04 MeV. The very good 
Xmin obtained, the good values for the q's, and this correct pion mass constitute convincing 
proof for the presence of chiral TPE loops in the long-range pp interaction. 

In order to investigate further the importance of %TPE, we move the boundary inwards to 
b = 1.4 fm. When only OPE is used as long-range force, it is possible to achieve a reasonable 
fit: at the cost of 31 BC parameters Xmm = 2026.2 is reached. We then add to OPE the 
potential TPE(l.o.). With 28 BC parameters, Xmm = 1984.7 is obtained. Compared to 
only OPE, this corresponds to a drop in Xmm °f 41.5 with 3 fewer parameters, a significant 
improvement. However, the fit is still not optimal [j34" 1. We next add also the potential 



TPE(s.o.). With fixed c x = -0.76/GeV, this gives with 23 BC parameters xLn = 1934.5, 
c 3 = -4.99(21)/GeV, and c 4 = +5.62(59)/GeV. This shows that OPE together with xTPE 
gives a very good NN force at least as far inwards as 1.4 fm. 

In conclusion, we have, for the first time, incorporated and studied chiral TPE in an 
energy-dependent PWA of the pp scattering data. The main result of this Letter is that we 
have shown the presence of chiral TPE loops in the long-range pp interaction. A significant 
improvement over using just OPE is seen. With OPE and %TPE, an excellent fit to the 
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database becomes possible, even somewhat better than the standard 1998 pp PWA. The 
chiral parameters agree with those found in ttN scattering. Especially important in obtaining 
the very good fit is the isoscalar central attraction from the C3-term, partly a "chiral van- 



der-Waals force" due to the axial polarizability of the nucleon ]35] . In all, our results provide 
a big success for chiral symmetry. A novel class of PWA has been established, with such 
a theoretically well-founded and model-independent chiral TPE potential included in all 
partial waves. 
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TABLES 



TABLE I. Coefficients of the subleading-order potential TPE(s.o.) of Eq. (^|), for the central, 
spin-spin, tensor, and spin-orbit terms, both isoscalar and isovector. We defined c*o = co/g\, 
di = CiM/g\ for i = 1, 3, 4, and c~04 = cq + 4c4. 
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TABLE II. Results for the PWA's with different long-range interactions. #BC is the number 
of BC parameters. 







b = 1.4 fm 
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#BC 
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Nijm78 
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